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Abstract 





HIBALL - A Conceptual Heavy Ion Beam Driven Fusion Reactor Study 





A preliminary concept for a heavy-ion beam driven inertial confi- 
nement fusion power plant is presented. The high repetition rate 
of the RF accelerator driver is utilized to serve four reactor 
chambers alternatingly. In the chambers a novel first-wall protec- 
tion scheme is used. At a target gain of 83 the total net electri- 
at output is 3.8 GW. The recirculating power fraction is below 

15 > e 


The main goal of the comprehensive HIBALL study (which is continu- 
ing) is to demonstrate the compatibility of the design of the 
driver, the target and the reactor chambers. Though preliminary, 
the present design is essentially self-consistent. Tentative cost 
estimates are given. The costs compare well with those found in 
similar studies on other types of fusion reactors. 


Zusammenfassung 





HIBALL - Eine Konzeptstudie für einen @urch Schwerionenstrahlen 








getriebenen Fusionsreaktor 





Es wird ein vorläufiges Konzept für ein Trägheitsfusionskraftwerk 
mit Schwerionenstrahltreiber vorgestellt. Die hohe Pulsfolge- 
frequenz des Treibers vom HF-Beschleuniger-Typ wird ausgenutzt, 
um vier Reaktorkammern abwechselnd zu bedienen. In den Kammern 
wird eine neuartige Technik zum Schutz der ersten Wand einge- 
setzt. Bei einem Target-Energiegewinn von 83 beträgt die elektri- 
sche Netto-Gesamtleistung 3.8 GW. Der rezirkulierende 
Leistungsanteil liegt unter 15%. 


Hauptziel der umfassenden HIBALL-Studie (die fortgesetzt wird) 
ist es, die Verträglichkeit der Entwürfe für Treiber, Target und 
Reaktorkammern nachzuweisen. Obwohl vorläufig, ist der vorliegen- 
de Entwurf im wesentlichen konsistent. Es werden vorläufige 
Kostenschätzungen angegeben. Die Kosten sind gut vergleichbar mit 
denen, die in ähnlichen Studien für andere Typen von Fusionsreak- 
toren gefunden wurden. 
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References are found after each subchapter (e.g., III.1) or chapter (e.g., I) 


























VI.1 General Requirements, Description, and Reference Parameters 





VI.1.1 Requirements 





The requirements imposed on the design of the cavity for HIBALL are 






determined by a number of sometimes conflicting considerations. A strict 
definition of the cavity is that part of the chamber inside the first material 


wall. For the purposes of this section the definition is expanded to include 





the entire chamber including reflectors, shields, and the top and bottom 






portions. The basic requirements are: (1) the cavity be able to dissipate the 
energy from the target in a coolant operating at temperatures consistent with 
material requirements and the necessity to produce steam suitable for oper- 
ation of a steam turbine driven generator, (2) the atmosphere in the cavity be 
such that the heavy icn beam can be focussed on the target, (3) the system be 
maintainable in a reasonable fashion, (4) the system be able to breed tritium 
at a breeding ratio greater than or equal to unity, (5) the tritium be re- 
coverable without undue accumulation of tritium in the breeding material, and 
(6) the dose rate outside the shielding be at an acceptable level. 

The two requirements which have the greatest impact are the need to 
dissipate the target energy and propagate the beam through the chamber to the 
target. The fusion yield of the HIBALL target is 400 MJ and the repetition 
rate is 5 Hz, yielding a fusion power of 2000 MW per chamber. The time de- 
pendent neutronics studies reported in section VI.3.4 indicate that because of 
endothermic reactions in the target the actual energy available is slightly 
less, i.e. 396 MJ/shot. Of these 396 MJ, 285 MJ is carried by neutrons of 


average energy of 12 MeV, 90 MJ is contained in X-rays emitted from the tar- 


get, and 21 MJ by the target debris, i.e. He, unburned D and T, Pb, and Li. A 


very small amount appears as high energy gamma radiation. This means that al- 
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though the greatest part of the energy is carried by the neutrons and there- 
fore is volumetrically deposited in the blanket materials, a significant 
amount, viz. 111 MJ, is in a form that would be deposited on or very near the 
surface of any material surrounding the target. 


VI.1.2 INPORT Concept 





A persistant technica! problem in the Inertial Confinement Fusion (ICF) 
field has been the protection of the first load bearing walls from target 
debris, x-rays, and neutrons. Various schemes have been proposed in the past: 
swirling liquid metal pools(1). wetted walls[2,3,4,5) magnetic protection(6), 
gaseous protection(7), dry wall ablative shields (8-11), and free falling 
Sheets of liquid metals(12-15), None of these schemes have been completely 
satisfactory, but the liquid metal protection scheme, HYLIFE‘ 15) | seems to be 
the best developed thus far. 

One disadvantage with the HYLIFE scheme is the disassembly of the liquid 
metal columns after each shot and the need to reestablish the stream before 
the next target can be injected. This results in small AT, high recirculation 
rates, and large pumping powers, especially for the Pb-Li alloys. Ideally, 
one would like to develop a scheme which would slow down the flow of liquid 
metal so that it can absorb the energy of several shots before exiting the 
reactor. 

Such a scheme was developed in late 1979 for the HIBALL project. The 
basis of this new design is the use of woven SiC tubes which are flexible, 
sufficiently strong, compatible with Pb-Li alloys used in HIBALL, and porous 
enough to allow liquid to cover the surface while the bulk of the fluid flows 
down the center of the tube (Fig. VI.1-1). This idea is called the INPORT 
concept, standing for the Inhibited Flow - Porous Tube Concept. The film 


thickness of roughly 1 m is sufficient to absorb the energy from x-rays and 
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Figure VI.1-1 
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target debris while several banks of tubes provide enough Pb-Li alloy to con- 


siderably moderate the neutron flux and reduce the total damage and damage 


rate in the first wall (Fig. VI.1-2). 


A schematic of how the porous tube structure fits into the HIBALL chamber 


design is given in Fig. VI.1-3 and the overall chamber design is shown in Fig. 


VI.1-4. A few of the major operating parameters are given below: 


Tube radius, cm (Number) - inner rows 


- Outer rows 


Inlet/outlet temp. of Pbg3Li)7 


Max. flow velocity (front tube) 

Target yield x frequency of shots 
Energy flux to first Pb-Li surface 
Condensation time to 10”? torr (at RT) 
dpa in HT-9 FW per FPY w/o INPORT units 
dpa in HT-9 FW per FPY with INPORT units 


1.5/ (1230) 
5/(3070) 
330/500°C 

5 m/s 

400 M) x 5 Hz 
35 J/cm 

0.15 sec 

21 

2.7 


It is clear that the INPORT concept enhances the promise of ICF reactors and 


it should remedy one of the more serious drawacks of the HYLIFE concept, the 


disassembly after each shot. The rest of this section will discuss, in more 


detail, how the INPORT concept fits into the HIBALL desian. 
VI.1.3 Mechanical Description of Cavity 





The cavity is an upright cylinder with internal dimensions of 11.5 m 


height on axis and 10 m in diameter, Fig. VI.1-4, 


From the center and outward 


in all directions, the cavity is characterized by three distinct zones, the 


blanket, reflector and shield. 


In the vertical sides of the cavity, the 


blanket consists of a 2 m thick zone of SiC tubes through which Liy7Pbg3 is 
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Fig. VI-1-2 Effect of INPORT tubes on Damage in HT-9 first wall. 
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circulated, Fig. VI.1-3. It is followed hy a reflector zone made of HT-9 
which is 40 cm thick and finally by the shield, which is primarily concrete 
and is 3.5 m thick. The top of the cavity has wedge shaped blanket modules, 
50 cm thick which are also filled with Li,7Pbg3. They too are followed by 40 
cm of HT-9 reflector and a 3.5 m concrete shield. The bottom of the cavity 
has a pool of Li17Pbg3 one meter thick, followed by a 40 cm thick reflector 
which has holes in it to allow the coolant to drain out. A catch basin then 
directs the flow to an outlet tube through which it is pumped to the steam 
generator. The rest of the material that follows is a concrete shield. 

The top of the cavity has a unique design. Although the blanket modules 
are designed to stay stationary, the reflector and shield is capable of being 
rotated about the central axis. This is needed to provide access to the 
blanket modules for replacement and will be discussed further in Chapter IX. 
The upper blanket moduies are locked to the reflector by means of studs which 
fit in milled slots. During reactor operation when the modules are filled 
with Li,7Pbg3, the stud latches are actuated and the modules are locked to the 
reflector. However, during servicing when the modules are drained of all the 
breeding material, the stud latches are de-activated. At this time the 
blanket modules are only supported on the outer periphery and in the center of 
the cavity, where they are attached to a central hub. This makes it possible 
to rotate the reflector and shield while the blanket modules remain station- 
ary. 

The cavity is sealed to the outside by a liner which is welded to the 
reflector at the upper end of the cavity. A circumferential weld (or seal) 
hetween the liner on the stationary part of the top shield and the rotatable 


part maintains the cavity atmosphere during operation. 
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At the iunction betveen the top blanket modules and the SiC tubes, there 
are 30 apertures, 65 cm high and 1.2 m vide. These are the vacuum system 
ports. The vacuum ducts which are concealed from direct line of sight of 
neutrons lead to pumping stations located in the upper corner of the reactor 
cavity. These pumps are used to evacuate the cavity prior to operation and to 
pump the non-condensable gases, such as the hydrogen species and the helium. 
The pumping system is discussed further in section VI.5. 

There are 20 beam ports in each cavity. Each two beams come in at + 16° 
to the horizontal, spaced at 36° circumferentially. The blanket is divided in 
such a way that every third module has two beam ports built in. At the first 
surface (5 m radius) the beam ports are 23 cm wide and 74 cm high. Distri- 
bution manifolds surround the beam tubes such that the SiC tubes are attached 
to them on the top and bottom. Each beam transport line is attached to the 


cavity at the point of entry and constitutes the vacuum interface between the 


cavity and the beam handling system. 


VI.1.4 Chamber Parameters 





Table VI.1-1 presents reference parameters which have been used as a 
basis for the calculations in the following sections of this chapter. These 
values were chosen on the basis of: (1) defining the nature and role of the 
system, (2) consistency with known materials characteristics of various parts 


of the system, or (3) as a result of previous calculations. A more detailed 


parameter set is given in Appendix A. 
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Table VI.1-1 





Cavity and INPORT Region 





Neutron energy/shot (MJ) 
X-ray energy/shot (MJ) 
Debris energy/shot (MJ) 
Gamma energy/shot (MJ) 
Cavity shape 
Cavity diameter to vacuum wall (m) 
Cavity height at center (m) 
Coolant and breeding material 
INPORT tube structural material 
INPORT region support structure 
First surface radius (m) 
INPORT tube region packing fraction 
INPORT tube length (m) 
INPORT tube diameter (cm) 
First two rows 
Remainder 
Number of first row tubes 
Number of remaining tubes 
Number of Beam penetrations 
Total area of beam penetrations at the 
first surface (m2) 
Pb atom density (x 10710 atom/cm 
just before shot 
Noncondensible atom density at 500°C 
(x 10-10 atoms/cm?) 
Chamber top thickness (m) 
Coolant volume fraction 
Bottom region thickness (m) 
Coolant volume fraction 
Vacuum Wall (first material wall) 
Structural material 
Inside diameter (m) 
Thickness (m) 


3) 





— Y — 


285 

90 

21 

< 1 
cylindrical 
14 

11.5 

Pbg3L1 37 (natural) 
SiC 

HT-9 

5 

0.33 

10 


0.01 





— — — — 
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Table VI.1-1 (continued) 


— D — 





Reflector 
Structural material 
Cool ant 
Coolant volume fraction 
Inside diameter (m) 
Thickness (m) 

Shield 
Structural Material 
Cool ant 
Coolant volume fraction 
Inside diameter (m) 
Thickness (m) 

Pbgaliız Cocl ant 
Inlet temperature (°C) 
Outlet tempe-ature (°C) 
Pressure (MPa) 














HT-9 
Pbg3Lil7 
0.9 
14.02 
0.4 


Concrete (unreinforced) 
0.05 

14.82 

3.5 


330 
500 
0.2 
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VI.2 Impulse Response of INPORT Tubes 





VI.2.1 Impulse Characterization 





When the PbLi is vaporized by the target generated X-rays, a reactive 
impulse imparts radial acceleration to the first two rows of the tube bank. 
One can conservatively estimate the bulk kinetic energy of the vaporized fluid 


to be equal to the thermal energy of the gas. 


(3/2) (M..../ KpTgas © (1/2) Mg, y (VI.2-1) 


gas Mion) B gas < 





Here the left-hand-side represents the thermal energy of the gas which is 





equal to the deposited X-ray energy less the energies of ionization and 





vaporization. The bulk velocity of the gas is denoted by V while Mion 





represents the average mass of ions in the gas. For the case of 87.6 MJ of X- 






ray energy and 13.3 kg of vaporized PbLi, Tgas = 1.26 eV. The reactive 





impulse equals MaasV- 















1/2 





AI = M (3k,T 


gas gas! ion) 








(VI.2-2) 
= 1.89 x 109 dyne-sec. 









This can also be expressed as an impulsive pressure per unit surface area. 








AP = 600 dyne-sec/cm° (VI.2-3) 








The preceding estimate is conservative since it is assumed that all vaporized 
atoms are initially moving towards the center of the cavity with their thermal 


velocity. This will clearly overestimate the reactive impulse. 
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V1.2.2 Stiffened INPORT Tubes 





V1.2.2.1 General Characterization 





Both stiffened and unstiffened silicon carbide tubes were considered for 
the INPORT units. In the former case the woven tubes are rigidized with a SiC 
binder or matrix but with a porous wall. Such members resist transverse 
loading by flexural action, not unlike structural beams. A general formu- 
lation is developed for stresses, deflections and frequencies. For uniform 
impulsive pressure, design data is presented which identifies the basic 
mechanical characteristics of these elements. 


VI.2.2.2 Dynamic Analysis 




















A solution of the form 





The principle equation of motion governing the transverse vibration of 


the tube is given by 


EI ə” y(x,t) /ax’ + m ə” y(x,t) /ət = p(x,t) (V1.2-4) 


where E = elastic modulus 


I cross section moment of inertia 


M = mass per unit length 


p(x,t) = dynamic pressure (force per unit length) 
t = time 
y(x,t) = transverse displacement. 
The coordinate system and loading are shown on a typical span of tubing in 


Fig. VI.2-1. 


To determine the natural frequencies of the tube, the external load is 


set equal to zero, resulting in the differential equation for free vibration. 





p(x,t) 





































ATT 














> + 
> 2 
Ax 
ຈ. 
> 
Fig. VI.2-1 


























mö 77 


y(x,t) = e(x) Y(t) (V1.2-5) 


is substituted into Eq. VI.2-4, leading to 


EI den) 5 vı, Cu (V1.2-6) 


mo(x) dx dt 





With x and t being independent variables, the preceding equation can be satis- 
fied only if both sides are set equal to a constant, producing two ordinary 
differential equations. Choosing u” for this constant, these equations will 


have solutions of the form 


Y(t) = ໃງ sin wt + Co COS ut (VI.2-7) 
(x) = A sin ax + B cos ax + C sinh ax + D cosh ax (VI.2-8) 
where 
aca rt > 


The tubes would be continuous through any intermediate support, if pro- 
vided, with the deformed shape for any span characterized by zero slope change 


and zero deflection at each end. Thus, the boundary conditions are 


y(o,t) = 7 (2,t) = 0 
(VI.2-9) 


aylo,t)/ax = ay(£,t)/ax = O . 








v1.2.5 —#— 
With these and Eq. VI.2-8 the following transcendental equation is obtained. 
605 81 cosh ag = 1 (VI.2-10) 


This equation generates an infinite number of discrete eigenvalues, 


an(n=1,2,3,...), which are obtained numerically and are related to the natural 


frequencies w, by 


w = aztel/m)1/7 ຄ = 1,2,3, ... (VI.2-11) 


Corresponding to each natural frequency is a characteristic shape or 


eigenfunction 


ໃ (x) = C„L(A,/B )(sirh a X - sin a x) + cosh a x - cos a,x] (VI.2-12) 


where 


A/B, = (cos a t - cosh a,£)/(sinh a £ - sin ant) 


and C, is an arbitrary constant. 


After the mode shapes and frequencies have been determined, a modal- 
superposition method of analysis can be used to solve for the dynamic response 
of the forced vibration. Thus, any displacement can be expressed by Super- 


imposing appropriate amplitudes of the vibration mode shapes for the 


structure. This can be written as 


y(x,t) = 2 ໃ (x) Y(t) (VI.2-13) 





| 
| 
| 
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Using an energy solution, Lagrange's equation is applied to the system 


d[əK/əY (t)1/dt + 3U/aY, (t) = an, /əY (t) (VI.2-14) 
where K = total kinetic energy 
U = total flexural strain energy 


Q. = external work 
and (*) denotes differentiation with respect to time, t. The kinetic energy 


for the entire system is given by 


>< 
H 


(1/2) m 139(x.y)” dx 
(1/2) m [EL] e (x) Y (t)J2 dx 
n 


(V1.2-15) 


= ç v2 £ 2 
1 /2 Y-(t d 
(1/2) m l nt) Jo 0) dx 


vhere m 1s considered to be constant over the span and orthogonality of the 
mode shapes has been utilized. 


The work done by external forces is 


p(x,t) y(x,t) dx 


: b (V1.2-16) 
= Je Pt) CE a(x) V(t] dx . 


Considering p(x,t) = P(x) f(t), where P(x) is the spatial load distribution 


and f(t) is the time function, the external work becomes 


ດູ = f(t) J Y (t) fo P(x) e (x) dx . 
n 
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VI .2-7 


Substituting into Eq. VI.2-14 yields 


a Y (t) JÉ 9ğ(x) dx = auzar (t) = F(t) Jo P(x) 8,00 dx 


or (VI.2-17) 


Y (t) = ağ Y (t) = f(t) (5 P(x) @,(x) ex/ñ Jó +C (x) dx . 


Eliminating the time dependency, the modal static deflection can be 


defined as follows. 










R à F. kd 
Y, i“ le P(x) &, (x) dx/m w be @ (x) dx . (VI.2-18) 
Yo determine the dynamic response, the static response is multiplied by the | | 
dynamic load factor (DLF). More descriptively this could be identified as a 


dynamic load function, representing the time function of the loading for an 


Sk — - — ——. 


equivalent single degree of freedom system. Additional details on this can be 


found in Reference (1). The total deflection at any point is given by 



















y(x,t) = ) & (x) Və (DLF) , ` (VI.2-19) 
n 









The dynamic bending moment can be determined from 











M(x,t) = EI 22y(x,t)/ə2 = EI EY, st do (x)/dé DLF, . (VI.2-20) 
n 









The flexural stess in the tube is 








o(x,t) = M(x,t) R /I = ER, ; k d° + (x)/dx” (DLF) . (V1.2-21) 
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Here R, represents the outer radius and the 


at the outer surface. 


v1.2.2.3 Quantitative Results 





The preceding work identifies principle 
response. Complete details were incorporate 
determination of numerical values. 

The fundamental frequency curves are sh 
for various radii and thicknesses in Figs. V 
the mass of the PbLi. It can be seen that 1 
frequencies are very low and relatively cons 
cm, the fundamental frequencies rise very qı 
results indicate that the effective length | 
siderably less than the full length of 10 m. 
results based upon a dynamic load factor of 
resonance is not considered. As the tube {| 
e.g., 1 mm x 1.5 cm, the transverse displacı 

effects can be expected for flexural stress: 
dicted for an INPORT unit with a radius, th 
and 10 m is unacceptably high. 

Natural frequencies can be raised subs 
while deflections and stresses are reduced 
nism is employed. Figure VI.2-7 is a conce 
support which sustains the radial impulsive 
cumferential tensile stresses. Such a desi 
transmit thrust to the outer wall. Small r 
sketch, would generally be important to red 





stress is associated with points 


equations for stiffened tube 
d into a computer program for the 


own as functions of tube length 
1.2-2, 2-3 and 2-4. These include 
or lengths greater than 800 cm, 
tant. For lengths less than 200 
ickly. To avoid resonance, these 
between supports) should be con- 

, Figure VI.2-5 shows deflection 
unity for a single impulse, i.e., 
nickness and inner radius decrease, 
ments rise very rapidly. Similar 
as. In Fig. VI.2-6, the value pre- 


ickness and length of 1.5 cm, 1 mm, 


tantially above the repetition rate 

if an intermediate support mecha- 
ptual representation for a ring 
forces by self-equilibrated cir- 
gn precludes the necessity to 

adial dimensions, as implied by the 


uce thermal stresses caused by 
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Fig. VI.2-2 


FUNDAMENTAL FREQUENCY 


VS. TUBE LENGTH 
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Fig. VI.2-3 


FUNDAMENTAL FREQUENCY 
VS. TUBE LENGTH 





MATERIAL : STIFFENED SiC 
COOLANT : PbLi 
TUBE THICKNESS : 2.0 mm 
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Fig. VI.2-4 


FUNDAMENTAL FREQUENCY 
VS. TUBE LENGTH 
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Fig. V1.2-5 


MAXIMUM DYNAMIC DEFLECTION 
VS. TUBE RADIUS 
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Fig. VI.2-6 


MAXIMUM DYNAMIC FLEXURAL STRESS 
VS. TUBE RADIUS 





MATERIAL: STIFFENED SiC 
COOLANT : PbLi 

TUBE LENGTH: 1000 cm 
DAMPING: 2% CRITICAL 
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radial temperature gradients. Active protective cooling would still be es- 
sential. 

A combination of reasons led the decision to use unstiffened rather than 
stiffened SiC tubes for the INPORT units. In the latter case, the geometry is 
that of a somewhat “elicate structure requiring extra care in fabrication and 
handling. The analysis has shown that resonance, large deflections and 
stresses are potential problems without supplementary structural support which 
of course is another addition to complexity. The stiffened tubes also involve 
untested manufacturing techniques and unknown additional expenses as well. By 
comparison, unstiffened tubes do not appear to have such disadvantages while 
possessing inherent desirable features. Their mechanical response character- 
istics are discussed in the work which follows. 


VI.2.3 Mechanical Response of Unstiffened INPORT Tubes 





VI.2.3.1 General Characteristics 





The term "unstiffened" used here denotes a tube woven from silicon car- 
bide fibers alone. The wall is pliable and porous, essentially cloth-like. A 
small percentage of fibers will run axially in the wall to enhance the tensile 
load carrying capacity. 

As structural members, such units resist lateral loading by tensile force 
action alone, i.e., having negligible shear and bending resistance. For the 
dynamic response, dissipation is accounted for by including viscous damping in 


the system. In the work which follows, two support modes are considered sepa- 





The first involves suspension from the upper end with the lower free 





rately. 


to move while the second has the bottom end constrained in addition. 


VI.2.3.2 Equilibrium Equations and Vibration Modes-One Point Suspension 





The notation and coordinate system for the equilibrium equations is shown 


in Fig. VI.2-8. On the element, T represents the tension while p and g denote 








INPORT TUBE EQUILIBRIUM NOTATION 
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Fig. VI.2-8 
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the mass density per unit length and the gravitational constant, respectively. 
At a generic location x, the lateral displacement component is taken to be y. 
For response characterized by small slope changes, equilibrium in the axial 
direction requires the tension gradient to be equal to the weight per unit 


length, 
dT/dx = pg 
ad (VI.2-22) 


T = pgx 


from the condition that the tension is zero at the lower end. In the trans- 


verse direction, the equilibrium statement can be reduced to 
T 30/ax + 099 = p 3”y/ət” (VI.2-23) 
in which t denotes time. The two preceding equations can be combined into 
X ə”y/ əx” + 9y/öx = ə”y/gət” ə (VI.2-24) 


Characteristic solutions may be developed by taking the displacement as a 


product of functions in the variables: X(x) T(t). Thus 
T/gT = xX"/X + X'/X = - 9% . (V1.2-25) 
It follows that the solution for the time-dependent component 1s 


T = A cos Agt + B sin avgt (V1.2-26) 
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where 1, A and B are constants. The equation for the spatial component is of 


N | 1/2, _ 
the form ກກ "າ q... For a single impulse the initial conditions are 





5 DEFLECTION 
y(x,0) = 0 ay(x,0) /at = f(x) (V1.2-29) i 3 cm TUBE , 


With these and orthogonality of the solutions, 


r-2, The roots are 


(x y")" + (ax + bx y = 0 


i 
SUE 


0.1877 Hz 0.4309 Hz 0.6756 Hz 


2,21/2 = 2,4048, 5.5201, 8.6537, etc. 





the solution of which can be expressed in terms of Bessel Functions of the 


first and second kind, J, and Y,. In general the parameters and solution are 








However, the natural frequencies can be defined as . fŠ 1/2 1/2 1 າ2 1/2 
„su Ñ. So f(x) ປງ (A! )dx/g An Je (217% )dx (VI.2-30) 


m 


a = (1-r)/2 Y = (2-r+s) /2 


=" agl”? 


| In this work, f(x) is taken as a uniform initial velocity Yo- Then after | | 


1/2 | 


/(2-r+s) evaluating the integrals in Eq. VI.2-30, the result is | 


k = 2la] 





and thus, for example, the fundamental is 





1/2 J 











2 1/2 
v = [.1-ກໍ່ - 4] '"/(2-r+s) . $ ne { 1/2 .2- 
1/2 B Ya An (g2) 22 18 ) (VI-2-31) 


(rad/sec) Ë 





ພາ = 1.2024 (9/2) 















































| 
y = x [CJ (kx7) + DY, (kx')J * ə lution is complete for the initial impulse problem. Numer- | 
I š From this it can be seen that system frequencies are determined essentially by NO. FREQUENCY (Hz) ee | ; 
ical results are obtained from a computer program based upon the preceding vən, 
w... s ww A A idin | 0.1877 analysis. Calculations include the weight of the PbLi, the SiC and the outer + | 
2 0.4309 surface Film. | 
pei b = 0 a = x s = 0 .” ະທ ee 3 0.6756 The displacement profile following an impulse at rest is initially very 
a= 0 Y = 1/2 k = 2 və . 4 09205 uniform. A typical result is shown in Fig. VI.2-10 at a time of 0.20 sec. | 
ກາບກ SEEN ຍາກຫ າແຕ RUN: ກະແທກ. ກກກ. ww. ` This is the instant at which a second impulse would occur in operation. How- | 
Thus Shapes corresponding to the first three natural frequencies are also shown. 5 1.1656 ever if the excitation is limited to a single shot, and motion allowed to de- 
V1.2.3.3 Forced Response - One Point Suspension 6 1.4519 velop for larger values of time, the shapes become somewhat complex as shown | 
F CJ, (2axl/2) ë DY (2xx1/2) (VI.2-27) The general form for the transverse displacement can be obtained using 7 | 7042 qualitatively in Fig. VI.2-11. The quantitative time histories of Figs. 
ERECTED ຕ ram ຫະກາ. R 19566 V1.2-12 and VI.2-13 correspond to this case also. The maximum midpcint and — > 
Since Y, becomes arbitrarily large as x tends to zero, it is necessary to 9 .— tip deflections are approximately 16 and 55 cm, respectively. It was felt 





delete this component, or simply set D equal to zero. Satisfaction of the 











O 


2.4613 














k. ` 
SAA 7:00 





n..” ””- 











sáb | 








































































































V1.2-23 Ë “.. : — 46 — ) “”— 
-22  -ຈ-- | Fig. VI.2-11 1.8 | VI.2-25 v1.2-26 i 
DEFLECTED SHAPES OF 3cm TUBES | | 
2-10 AT VARIOUS TIMES AFTER UNIFORM IMPULSE — vs values, these displacements were excessive, particularly on a 
single shot basis. 
£ VI.2.3.4 Analytical Outline - Two Point Suspension 
PROFILE OF 2.5 $ 5 S “5 S IO 5 E ແ z 5 “es In order to control deflections, it has been proposed to constrain the 
ມ Prr: >. PR INPORT TUBE MECHANICAL RESPONSE ມ io Lx S : ກະບ 
AT t ກ 0.20 sec | | ə a ລ ç; 5 K... bottom as well as the top of the INPORT tubes. The background analysis for 
$ IMPULSE : 600 dyne-sec/cm® 2 2 645 ມ this case is briefly outlined in the work which follows. 
ໂວ TUBE : 3 ເຕ x IO m O | oz £ 4 | ມ 
| | = 15 x ve DEFLECTION:15.7 cm en si EOS a E The principle equation of motion is 
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i.e., b, (o) and b,(0) are the Fourier coefficients of f(x), g(x). With the 


preceding, the solution tc Eq. VI.2-32 can be expressed as 


y(x,t) = e bt ) Cf, cos ພາ + (9, + bf )/(sin wt) /w, Jsin nnx/£ . (VI.2-35) 





For a single impulse, with no initial displacement 






















ay(x,0)/at = g(x) = I(x)/p 
where I(x) is the impulse per unit length. Thus Eq. VI-2-35 becomes 


sino t 
y(x,t) = (28% for) E [J 1(u) sin(mmu/2)du] — sin X. (VI-2-36) 








If a more general force is imposed, F(x,t) per unit length, 


2e-bt w (t-r) 


£ pt . NTU ... nx: n 
y(x,t) = 7 ) le je F(u,t) sin — sin — sin — drdu . 








/ [ . 2 .3 .5 Numerical Results - Two Point Suspension 





The response of the INPORT tubes was first determined for a single im 
pulse, uniformly distributed along the length. Figures VI.2-14 - VI.2-17 show 


the time histories of the midspan displacement for various levels of damping. 





The amount of damping is expressed as a fraction of the so-called critical 





value, i.e., the damping coefficient which would produce non-oscil latory 





response. The peak values shown for a single pulse appear acceptable. 






The program was extended to determine INPORT tubes response under repet- 









Figure VI.2-18 corresponds to "start up", with the 








itive impulses at 5 Hz. 










initial velocity and displacement zero. At the end of 0.20 sec the velocity 


Fig. VI.2-14 
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and displäcement functions are evaluated and used as starting conditions with 
another uniform impulsive pressure. The process is repeited and carried out 
for fifty shots. These curves are for the first fifteen firings and show a 
small amount of overshoot with eventual stabilization about steady state mean 
values. On a relative basis the steady state mean amplitudes are moderate, 
e.g., less than 2% of the length. The limiting case of zero damping is shown 
for comparison in Fig. VI.2-19 along with the 20% value. The latter is recom- 
mended as a conservative estimate for such units. While theoretical and 
experimental values for metallic reactor components are much smaller single 
digit percentages, here the tube has a significant natural damping capacity 
from the woven soft construction and the ever present PbLi in liquid form. A 
level of 20% is a representative design figure for long span electrical con- 
ductors and transmission lines and is typical of experimental measurements 
(Ref 2). 

Finally the response was determined following shutdown. Typical results 
are shown in Fig. VI.2-20 in which the load sequence was stopped after fifty 
shots. While the displacement history is very smooth, characteristic of a 
Single mode, it should be noted that this result involves the contributions of 
the first eight modes in the series. Such results confirm expectations that 


no unusual characteristics develop for this scenario. 
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VI.3 Neutronics and Photonics 





VI.3.1 Introduction 





The neutronics and photonics analysis for HIBALL is presented in this 
section. A one-dimensional spherical geometry analysis was performed to 
determine the optimum blanket design which gives the highest possible energy 
multiplication with adequate breeding ratio. The blanket consisting of the 
INPORT tubes is also required to protect the ferritic steel structure from 
radiation damage. The biological shield is also designed. The actual 
cylindrical cavity is modelled in a detailed three-dimensional neutronics and 
photonics analysis. All cavity geometrical details are included in this 
analysis. A time-dependent neutronics analysis is also presented to account 
properly for the pulsed nature of the neutron source. This analysis yields 


time-dependent radiation damage rates which form an important input to micro- 





structure evolution calculations. Finally, the calculations of blanket and | 
shield radioactivity and afterheat are given. 


V1.3.2 One-Dimensional Time Integrated Studies 





VI.3.2.1 Introduction 





Neutronics and photonics calculations are required to determine important 
reactor parameters such as tritium breeding, nuclear heating, and radiation 
damage. In an inertial confinement fusion reactor, the DT fuel is heated and 
compressed to extremely high densities before it ignites. Neutron fuel inter- 
actions result in spectrum softening, neutron multiplication, and gamma pro- 
duction. In this section, a consistent coupled target-blanket neutronics and 
photonics study of HIBALL is given. The neutron and gamma source for blanket 
calculations was obtained from the target neutronics and photonics results 


given in section III.1. 








The HIBALL blanket is required to breed tritium and convert the kinetic 





energy of the fusion reaction into heat. Furthermore, the INPORT tubes are 


required to protect the vacuum wall and HT-9 structure in the reflector from 





radiation damage. One-dimensional steady state calculations are performed to 
determine the optimum blanket thickness which yields the largest possible 
energy multiplication with adequate tritium breeding ratio and first wall pro- 
tection. The thickness of biological shield required to reduce the biological 
dose to permissible levels is also determined. 


V1.3.2.2 Blanket and Shield Model 








The HIBALL fusion reactor power plant design incorporates four cylindri- 
cal cavities each having a radius of 7 m and a height of 10 m. The blanket 


region consists of an array of porous tubes (INPORT tubes) made of braided SiC 





through which the Li,7Pbg3 liquid metal eutectic flows. These tubes are 
placed inside the reactor cavity. Besides serving as the coolant and tritium 
breeder, the Li]7Pbg3 coolant is utilized for first wall protection. The 
tubes occupy 33% of the 2 m thick blanket region yielding an effective blanket 
thickness of 66 cm. The LiPb coolant occupies 98% of the tube volume with the 
SiC occupying the rene ining 2%. SiC and Li;7Pbg3 have nominal densities of 
3.17 and 9.4 g/cm, respectively. The first wall is made of ferritic steel 
(HT-9) which is an alloy composed of 85.25 wt% iron and 11.5 wt% chromium. 

The first wall has a thickness of 1 cm and a density of 7.3 g/cm. A 0.4 m 
thick reflector composed of 90 v/o ferritic steel structure and 10 v/o 
L117Pbg3 coolant is used. The reactor utilizes a 3.5 m thick concrete shield. 
The shield consists of 95% ordinary concrete and 5% H,0 coolant. A schematic 
of the blanket, first wall, reflector, and shield configuration for the HIBALL 


reactor is given in Fig. VI.3-1. 
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HIBALL reactor cavity model, 
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V1.3.2.3 Calculational Method and Nuclear Data 











The one-dimensional discrete ordinates code ANISN(1) was used to perform 
neutronics and photonics calculations yielding average time integrated re- 
sults. A P3-S4 approximation was used. Spherical geometry was used in the 
calculations and hence the results represent the conditions at the central 
plane of the reactor. The results give conservatively high damage rates and 
low tritium breeding ratios. The neutron and gamma spectra obtained from the 
target calculations are used to represent the source for the blanket neutron- 
ics and photonics calculations. The source is considered to be an isotropic 
point source at the center of the cavity. The calculations account for 
neutron spectrum softening, neutron multiplication, and gamma production in 
the target. The one-dimensional finite element code ONETRAN (2) implemented at 
Karlsruhe Nuclear Laboratory was also used to perform the calculations for the 
final blanket design and the results are compared to the results from ANISN. 

A coupled 25 neutron-21 gamma group cross section library was used. The li- 
brary consists of the RSIC DLC-41B/VITAMIN-C data library (3) and the DLC-60/ 
MACKLIB-IV response data library.(4) This data library is based on ENDF/B-IV. 
The nuclide densities used in the calculations are given in Table VI.3-1. The 
results presented here are based on a DT yield of 400 MJ and a repetition rate 
of 5 Hz yielding 7.1 x 1020 fusion neutrons per second. 

Neutron and gamma cross section requirements for the neutronic analysis 
of the HIBALL system are similar to those of fast fission reactor systems. 
Neutron data are required up to the fusion neutron energy of 14.1 MeV. Pre- 
sently available neutron cross section libraries like ENDF /B(5) and KEDAK (6) 
cover this energy range. Although data in these libraries show a lower level 
of confidence above a few MeV than in the keV and thermal energy range, and 


some improvement in data, e.g., energy distribution of secondary neutrons is 
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Table VI.3-1 
Nuclide Densities Used in Neutronics Analysis of HIBALL Fusion Reactor 


conceivable, the quality of the data as stored in ENDF/B-IV is considered, 












with exceptions given below, to be adequate at this stage of the study. 






One major parameter to be determined in the neutronic calculation is the 




























Constituent Nuclide Density 
Region Elements (nuclei/b cm) | tritium breeding ratio. Tritium is produced by the reactions 
Blanket Gli 0.00040 y 6 
[98 v/o Li, 7Pbaş ni 0.00505 odla — | 
+ 2 v/o SiC] Si 0.00095 
(.33 density factor) - 0.00095 The division of the tritium breeding between these two processes is essential- 
Pb 0.02661 














ly design dependent. Ströng doubts about the validity of the /Li(n,n'a)T 





















Ferritic Steel Fe .07171 cross section in ENDF/B-IV have been expressed, (7-10) Lowering the cross 
First Wall Cr .01039 section by 15-20% has been suggested. (8) Results of new differential measure- 
(1.0 density factor) Ni .00040 


00049 ments (Fig. VI.3-2) from different laboratories suggest reducing the ENDF/B-IV 


.00028 
.00042 
.00043 
.00078 
.00013 










cross section by about 10%. Some of these results are still of a preliminary 





nature. 







In the HIBALL blanket design, due to the presence of a large amount of 






lead, neutrons are rapidly slowed down below the threshold of the 7Li(nyn'a)T 












reaction. The contribution of this reaction to the tritium production is only 



























Reflector Fe 0.06454 
R 3 | ; i 
[90 v/o ferritic Cr 0.00925 about 2%. Thus ENDF/B-IV cross sections for this process can be used without 
Steel + 10 v/o Li; 7Pbg3] Ni 0.00036 | introducing significant error in the tritium breeding ratio. 
La 
(1.0 density factor) Mo 0.00044 The other reaction responsible for tritium breeding; 6Li(n,a)T 1s used in 
V 0.00025 11 
Si 0.00038 neutron data measurements as a standard. Recently Hale(11) has re-evaluated 
Mn 0.00038 | this cross section below 1 MeV for the ENDF/B-V standard file. He made a 
Ç 0.00070 thorough analysis of the ük compound nucleus system including charged parti- 
A 0.00011 
+; i i w i y R | 
TE 0.00004 cle reactions. His result (Fig. 11.3-3) gives a lower °Li(n,a)T cross section 
ük 0.00051 for the 252 keV resonance and for lower energies. For Gi, group cross 
0 


sections were generated using the ENDF/B-V Data and used to investigate .the 






effect on HIBALL overall tritium breeding ratio. 
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Figure VI.3-3 V1.3.2.4 Blanket Optimization 
ແ) 
I I | | pae A series of one-dimensional neutronics and photonics calculations was 
made in which the packing fraction for the INPORT tubes was varied with the 
2 ICO: bi - y | 1.0E +01 rrm TTTTTTİ | 
ລາ o * 7 blanket region thickness fixed at 2m. The effect of the tube packing 
N — 
LL ! = oa 2 . . ° ° . ° . p 
3 ə Jm ë > LI 6(N, ALPHA)T —— ENDF/BIZ fraction or the effective blanket thickness on radiation damage, tritium 
a — iz 2 i ) À ເກ 
breeding and nuclear heating was investigated. The design criterion used in 
ມ a —ENDF/BY . , , , ' 
a ə | + ເ ak the optimization was that the energy deposition should be as high as possible 
2 າ with a tritium breeding ratio > 1.1 in the approximate one-dimensional calcu- 
= À — = 
é y p = 5 lation. 
x 8 o > ທ Table VI.3-2 gives the effect of packing fraction on tritium production. 
=-= H—n— N ຈ Y ທ 
ro .. R 13 = o It is clear that tritium production in the blanket increases as the thickness 
< Z ul! u! < [ມ : ! * kə d | 
b e -<991 |). = > | increases. On the other hand, tritium production in the reflector region 
+, ul € I decreases as the effective blanket thickness increases. The net effect is 
x - © 5 = | | 8 ຈ້ | that the overall breeding ratio increases nearly linearly with the effective 
N | | 
ເວລາກ - w 6 š | blanket thickness. 
F > və ຈາ D Figure VI.3-4 gives the spatial distribution of tritium production in the 
ນ 
ຈ < D blanket and reflector for a packing fraction of 0.33. The contributions to 
* + © 3 
HO 8 tritium production from /Li(n,n'a) and 6Li(n,a) reactions are shown. The 
b he Jo dä contribution from Ni, which is a high energy reaction (threshold energy = 
E É 2.86 MeV), decreases sharply as one moves into the blanket away from the 
a > 
P າໃ ກ - source because of the increased softening of the spectrum. Table VI.3-3 gives 
\ > 
7 ENERGY. EV. a comparison between the results of ANISN and ONETRAN. The effect of using 
CREA Y “əli 
ດລ + m N - LOE-O1 ແທກ | L-I LI TM IL I I 11111 the ENDF/B-V data for the 6Li(m,a)T cross section is also given. Using the 
O O O O O 
1.OE+O3 1.0E+04 1.0E+05 1.0E +06 : 6Li(n,a)T cross section from ENDF/B-V results in a 2% reduction in the 











calculated tritium breeding ratio. 
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6, . i 
Li(n,a)T Cross Sections from ENDF/B IV and V Table VI.3-4 gives the neutron, gamma, and total energy deposition in the 





different regions for different values of the packing fraction. Notice that 
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Table VI.3-2 Effect of Tube Packing Fraction on Tritium Breeding 























Packing Fraction 0 .30 .33 -35 .40 

Blanket 0 1.094 1.177 1.226 1.331 
Reflector 0.043 | 0.022 0.018 0.016 0.012 
System Total 0.043 | 1.116 1.195 1.242 1.343 














Table VI.3-3 Comparison Between Results for Tritium Production 





Table VI.3-4 Effect of INPORT Tube Packing Fraction 





on Nuclear Energy Deposition (MeV/Fusion) 




















Tritium Production (Tritons/Fusion) 






























































Packing Fraction 0 ~ 30 .33 .35 .40 
Neutrons - -- 7.026 7.467 7.729 8.278 
Blanket Gamma — 6.663 6.687 6.697 6.705 
Total - - - 13.689 14.154 14.426 14.983 
Neutrons 0.407 0.020 0.016 0.014 0.009 
First Wall Gamma 1.141 0.203 0.180 0.163 0.125 
Total 1.548 0.228 0.196 0.177 0.134 
Neutrons 2.626 0.287 0.232 0.202 0.142 
Reflector Gamma 14.402 4.013 3.367 2.985 2.185 
| Total 17.028 | 4.300 | 3.599 | 3.187 | 2.327 
| | 
System Total 18.576 18.217 | 17.949 17.790 17.444 
Overall Energy 
1.326 1.300 1.281 1.269 1.245 


Multiplication 
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ANISN, SyP3s | ONETRAN, SgPz, | ONETRAN, SgPz> 
ENDF/B-IV ENDF/B-IV ENDF /B-V 
Region Element ÖL i(n,a)T 6Li(n,a)T SLi(n,a)T 
6} ; 1.1500 1.1558 1.1327 
Blanket ni 0.0270 0.0271 0.0271 
Total 1.1770 1.1829 1.1598 
First Wall Total 0.0 0.0 0.0 
6, ; 0.0180 0.0180 0.0175 
Reflector ni 0.0 0.0 0.0 
Total 0.0180 0.0180 0.0175 
Tritium Breeding Ratio 1.1950 1.2009 1.1773 
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more than 50% of energy deposited comes from gamma-ray heating because of 
neutron reactions in the HT-9 structure. The contribution from gamma heating 
decreases as the INPORT tube region thickness increases. As the blanket 


thickness increases, the total energy deposited in it increases with a larger 
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| ' 
~ fraction of this energy being deposited by neutrons. The larger energy depo- 


sition results from the larger nmber of neutron and gamma interactions taking 
+ place in the thicker blanket. The relative gamma contribution to the heating 
decreases as the thickness increases because a smaller portion of the blanket 
will be adjacent to the first wall where large gamma production occurs. The 


neutron and gamma heating in the first wall and reflector decrease as the 


blanket thickness increases because of the increased neutron attenuation in 


bititil 


the blanket. The total energy deposited in the system 1s Tound to decrease as 


l 


the blanket thickness increases. This results because of the decreased gamma 


| production in ferritic steel due to the lower neutron flux in the first wall 





and reflector. A neutron absorbed in ÖL: releases — 4.8 MeV while if it is 
absorbed in the structure - 7 MeV is released. The energy multiplication is 
also included in Table VI.3-4. This factor is defined as the total energy 


deposited in the system, including the energy deposited by X-rays and target 


a | 1 pagal 


debris at the first surface of the blanket, divided by the fusion reaction 
yield of 17.6 MeV. 


The spatial variation of the total power density in the blanket for a 
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packing fraction of 0.33 is given in Fig. VI.3-5. The power density decreases 
as one moves away from the source inside the blanket because of the increased 
neutron and gamma attenuation. The iarge gamma production in HT-9 causes the 


gamma heating to increase sharply as one approaches the first wall. This is 


Figure VI.3-4 Tritium production in the blanket and reflector. 


the reason for the existance of a minimum in the power density at ~ 30 cm from 


the first wall. 
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Figure VI.3-5 Spatial variatior 
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It is clear from the previous results that increasing the INPORT tube 


packing fraction results in increasing the tritium breeding ratio at the ex- 


pense of power production. Since energy production is the primary goal of a 


fusion power reactor, a packing fraction of 0.33 is chosen. This yields an 
energy multiplication of 1.28 and a breeding ratio of 1.195. A relatively 
high breeding ratio is chosen here because of expected reduction in tritium 


production when one includes the penetrations and various structural support 


materials in the detailed three-dimensional cylindrical geometry. The breed- 


ing ratio should also be greater than one to account for uncertainties in data 


and calculational methods. 


V1.3.2.5 Neutron and Gamma Spectra 








Using the neutron and gamma spectra from the target calculations to 





represent the source for the blanket neutronics and photonics calculations, 





the neutron and gamma spectra were calculated at different positions in the 





blanket, first wall and reflector. The neutron spectra are given in Fig. 





VI.3-6. In the innermost part of the blanket the spectrum has a pronounced 





peak at 14.1 MeV with the lower energy part of the spectrum resulting from 





neutron target interactions and neutron slowing down in the blanket itself. 





In the first wall and reflector the spectrum is considerably softened primari- 





ly because of slowing down in the blanket - reflector system. In Fig. VI.3-7 





the gamma spectra in the first wall and the center of the blanket are com- 





pared. The gamma photon density in the first wall is larger than that in the 





blanket because of the large gamma production in iron. 


VI.3.2.6 Radiation Damage to Structural Materials 





The effect of blanket thickness on the maximum number of atomic displace- 


ments (dpa) per full power year in the SiC INPORT tubes and the HT-9 first 


wall is given in Table VI.3-5. The average dpa rate in the SiC tubes 
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Table VI.3-5 Effect of Blanket Thickness on Maximum DPA Rate (dpa/FPY) 
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Effective Blanket Thickness 0 60 66 70 80 
(cm) 

SiC INPORT Tubes --- 116 118 119 121 

HT-9 First Wall 25.360 3.411 2.688 2.289 1.525 
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Figure VI.3-6 Neutron spectrum in the blanket, first wall, Figure VI.3-7 


and reflector. 


Gamma spectrum in the blanket and 
first wall. 
(The ordinate is the integral in 


(The ordinate is the inte each energy group.) 


gral in each energy 









increases slightly as the blanket thickness increases while the dpa rate in 
the ferritic steel first wall decreases as the bianket thickness increases. 
The reason is that the dpa cross section for iron peaks at 14.1 MeV and has a 
threshold energy of ~ 1 keV with the total dpa in iron decreasing as the 
spectrum becomes softer. On the other hand, the dpa cross section for 1 
peaks at ~ 3.5 MeV and has a lower threshold energy of ~ 0.2 keV resulting in 
an increase in the total dpa in C as the spectrum softens. 

The helium and hydrogen gas production rates in the HT-9 first wall are 
given in Table VI.3-6 for different effective blanket thicknesses. The aver- 
age gas production rate is found to decrease 85 the blanket thickness in- 
creases. The reason is that the hydrogen and helium production cross 
sections peak at 14.1 MeV. Furthermore, the reactions resulting in helium and 
hydrogen production are of a threshold nature with threshold energies in the 
low MeV energy range. The effect on gas production is more pronounced than 
the effect on dpa because of the lower dpa threshold energy. It is clear from 
the results that using INPORT tubes with an effective thickness of 66 cm re- 


duces displacement damage and gas production in the HT-9 first wall considera- 


bly, allowing it to last for the whole reactor lifetime. 
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Table VI.3-6 Effect of Blanket Thickness on Gas Production Rate 





in First Wall (appm/FPY) 











Effective Blanket Thickness 0 60 66 70 80 


(cm) 











Helium 229.1 0.659 0.364 0.245 0.090 


Hydrogen 





























Because of the large difference between the threshold energies for helium 
and dpa reactions, the helium to dpa ratio, which is an important parameter in 
microstructure calculations, decreases with spectrum softening. Figure VI.3-8 
Shows the effect of effective INPORT tube region thickness on dpa, helium pro- 
duction, and helium to dpa ratio in the HT-9 first wall. It is clear that the 
wall protection has a more pronounced effect on helium production. Conse- 
quently, the helium to dpa ratio decreases from a value of 9 for the unpro- 
tected wall to a value of 0.135 when an effective thickness of 66 cm is used. 
The helium to dpa ratio was also found to decrease as one moves inside the 
reflector. Figure VI.3-9 shows the spatial variation of dpa, helium pro- 
duction and helium to dpa ratio in the HIBALL reflector. The helium to dpa 
ratio decreases from 0.135 at the first wall to 0.019 at the back of the 
reflector. 

A comparison between the results for the HT-9 first wall radiation damage 
obtained using the ANISN coue and the results obtained using the ONETRAN code 
implemented at KfK, is given in Table VI.3-/. The differences obtained are 


attributed to the different methods used in the codes. ANISN treats the 
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Figure VI.3-8 Effect of INPORT tubes on Damage in HT-9 first wall. 
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Figure VI.3-9 Spatial variation of damage in HT-9 structure in HIBALL 


reflector. 


FIRST WALL 
REFLECTOR — 
L 1 1 1 
TOO TIO 720 730 7 








V1.3-22 — 


Table VI.3-7 Comparison Between First Wall Damage Results 





Of ANISN and ONETRAN 

















(S4P) 
dpa/FPY He appm/FPY H appm/FPY 
ANISN 2.688 0.364 1.380 
ONETRAN 2.700 0.319 1.230 














Table VI.3-8 Concrete Composition 











Element Weight % 
Calcium 24.33 
Silicon 20.24 
Oxygen 46.70 
Carbon 4.92 
Hydrogen 3.81 


























spatial variable in a finite difference form, while ONETRAN uses a finite 
element technique. 


VI.3.2.7 Shield Design 





Shielding of the reactor cavity does not pose special problems because of 
the simple geometrics involved as compared to shielding of the final focusing 
magnets discussed in section V.6. The reactor cavity shield is designed such 
that the biological dose rates outside the shield do not exceed ~ 5 mrem/hr 
during reactor operation. This design criterion has been set in order to 
allow hands-on maintenance of auxiliary components outside the reactor during 
operation. The shield for HIBALL consists of 95 v/o ordinary concrete (Type 3 
concrete from reference 12) and 5 v/o water coolant. Ordinary concrete was 
considered because of its relatively low cost and the absence of materials 
which can cause high activation. The concrete used has a density of 2.3 g/cm? 
and the composition shown in Table VI.3-8. 

A series of one-dimensional neutronics and photonics calculations was 
made to determine the flux attenuation as a function of shield thickness. 


Since a spherical geometry was used in the transport calculations, the results 

















correspond to the conditions at the reactor midplane where the dose is ex- 


pected to be greatest. Figure VI.3-10 shows the variation of neutron and 
gamma energy fluxes, at the shield outer surface, with shield thickness. It 
is clear that the neutron flux attenuates in the shield faster than does the 
gamma flux. This is due to the lack of high Z materials, which are more 
effective in attenuating gamma rays. 

The neutron and gamma fluxes at the outer surface of the shield were used 
to determine the corresponding biological dose values for different shield 
thicknesses. The tissue composition used in these calculations is given in 
Table VI.3-9. The variation of the anticipated dose rate with shield thick- 
ness is shown in Fig. VI.3-11. A shield thickness of 3.5 m results in a dose 
rate of 2.64 mrem/hr which is less than the design limit of 5 mrem/hr set in 
this work. Almost all of the contribution to the dose comes from gamma radi- 
ation. The shielding effectiveness can be improved by using high Z materials 
such as lead to attenuate the gammas. Boron which has a high thermal neutron 


absorption cross section with only a scft 0.48 MeV gamma being produced can 
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Figure VI.3-10 Variation of energy flux with shield thickness. 
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Figure VI.3-11 Variation of biological dose rate with shield 
thickness. 
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Table VI.3-9 Tissue Composition 























Composition Weight % 
Oxygen 76.2 
Carbon 11.1 
Hydrogen 10.1 
Nitrogen 2.6 





However, cost considerations become an important factor if 


A 3.5 m thick 


also be used. 


ordinary concrete is to be replaced by a more effective shield. 


concrete shield is used in HIBALL, resulting in an extremely low dose rate of 





2.64 mrem/hr outside the shield during reactor operation. 
Figures VI.3-12 and VI.3-13 show the neutron and gamma spectra at differ- 


ent positions in the shield at the reactor midplane. It is clear that signi- 


ficant neutron and gamma spectrum softening occurs in the first 0.5 m of the 


shield. At this point most of the neutrons are in the thermal group. The 


neutron flux attenuates as one moves deep in the shield with the spectrum 
shape remaining nearly the same. The gamma spectrum assumes a nearly fixed 
shape peaking at ~ 0.2 MeV at 0.5 m inside the shield. No further significant 
spectrum softening occurs as one moves deep in the shield. Again, the gamma 
flux does not attenuate as fast as does the neutron flux because of the gamma 
production following neutron capture and the lack of high Z materials. 

Figure VI.3-14 shows the spatial variation of power density in the re- 


actor midplane. The contributions from neutron and gamma heating are also 


shown. The peak power density in the shield is 0.045 W/cm®. The average 
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Neutron Spectrum in HIBALL Shield 
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Figure VI.3-12 Neutron s 


pectrum in HIBALL shield. 
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Figure VI.3-14 Spatial variation of power density at HIBALL midplane. 
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power density in the shield at the reactor midplane is found to be 1.83 x 1073 
W/cm. 


VI.3.2.8 Summary 





A one-dimensional spherical geometry neutronics and photonics analysis 
was performed to design the blanket and shield for HIBALL. The results show 
that by increasing the INPORT tube packing fraction, the tritium breeding 
ratio increases at the expense of energy multiplication. An effective tube 
region thickness of 66 cm was chosen resulting in an overall energy multipli- 
cation of 1.28 and a breeding ratio of 1.195. The INPORT tubes reduce radi- 
ation damage in the HT-9 first wall considerably allowing it to last for the 
whole reactor lifetime. A 3.5 m thick concrete shield was chosen resulting in 


an extremely low dose rate of 2.64 mrem/hr outside the shield during reactor 


operation. 


VI.3.3 Three-Dimensional Time Integrated Studies 





VI.3.3.1 Introduction 





The one-dimensional spherical geometry calculations presented in section 
VI.3.2 are useful for survey studies and predicting conditions at the reactor 
midplane. However, the one-dimensional calculations are not capable of 
adequately modeling the HIBALL reactor in which a point source is at the the 
center of the reactor cavity which has a cylindrical shape with a conically 
Shaped blanket at the top and a Li] 7Pbg3 pool at the bottom. Furthermore, the 
HIBALL blanket/shield system is required to accommodate large penetrations for 
vacuum pumping „nich cannot be modeled in a one-dimensional analysis. These 
geometrical effects are expected to have an impact on the system's overall 
tritium breeding and energy multiplication. In this section, a three- 


dimensional neutronics and photonics analysis is presented which adequately 


models the HIBALL reactor cavity. 
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The neutronics and photonics problem for an inertial confinement fusion 
reactor is basically a time-dependent problem. However, the concern here is 
only with time-integrated quantities such as the tritium breeding ratio and 
the total nuclear heating which can be obtained using steady state transport 
calculations. In this work, neutron multiplication, spectrum softening, and 
gamma production in the target is accounted for by using the results of the 
target neutronics and photonics calculations presented in section 111.1 as a 
source for the blanket and shield calculations. 


V1.3.3.2 Reactor Geometrical Model and Method of Calculation 





Figure VI.3-15 gives the blanket and shield configuration for HIBALL. 
The INPORT tube region has a thickness of 2 m and a height of 10 m. The 
INPORT tubes have a packing fraction of 0.33 with Lij7Pbg3 occupying 98 v/o of 
the tubes and the remaining 2 v/o occupied by SiC. The tube support structure 
is made of HT-9 and occupies 0.7 v/o of the tube region. A one meter deep 
Li17Pbg3 pool exists at the bottom of the reactor cavity. The upper blanket 
consists of 30 segments. The design of an upper blanket segment is shown in 
Fig. VI.3-16. The upper blanket consists of 97 v/o Li]7Pbg3» 1 v/o HT-9, and 
2 v/o SiC and is 50 cm thick. The Lij7Pbg3 in the region connecting the top 
blanket with the INPORT tubes helps protecting the HT-9 structure between the 
vacuum ducts. Thirty vacuum pumps are used to maintain the cavity pressure at 
1074 Torr. Each vacuum duct is 1 m wide and 0.6 m high. The first wall is 1 
cm thick and is made of HT-9. The side wall is 12 m high. The top liner is / 
and 6 m above the midplane at reactor centerline and side wall, respectively. 
A 40 cm thick reflector consisting of 90 v/o HT-9 and 10 v/o Li, 7Pbg3 is used. 
A 40 cm thick splash plate is used at the bottom of the reactor and is 
referred to as the bottom reflector in the following analysis. A 3.5 m thick 


biological shield made of 95 v/o ordinary concrete and 5 v/o H,0 coolant 
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